Extending electron orbital precession to the molecular case: Can orbital alignment be 

used to observe wavepacket dynamics? 
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The complexity of ultrafast molecular photoionization presents an obstacle to the modelling of 
pump-probe experiments. Here, a simple optimized model of atomic rubidium is combined with 
a molecular dynamics model to predict quantitatively the results of a pump-probe experiment in 
which long range rubidium dimers are first excited, then ionized after a variable delay. The method 
is illustrated by the outline of two proposed feasible experiments and the calculation of their out- 
comes. Both of these proposals use Feshbach 87 Rb2 molecules. We show that long-range molecular 
pump-probe experiments should observe spin-orbit precession given a suitable pump-pulse, and 
that the associated high-frequency beat signal in the ionization probability decays after a few tens 
of picoseconds. If the molecule was to be excited to only a single fine structure state state, then 
a low-frequency oscillation in the internuclear separation would be detectable through the time- 
dependent ionization cross section, giving a mechanism that would enable observation of coherent 
vibrational motion in this molecule. 

PACS numbers: 34.50.Gb, 33.80.Eh, 32.80.Fb 



I. INTRODUCTION 

The use of sub-picosecond laser pulses have enabled 
the observation of rapid processes in the time domain 
through the use of a pump-probe experiment. A typical 
experiment involves applying a "pump" pulse to a sam- 
ple, for example a gas of atoms or molecules, to initiate a 
time dependent process. A second "probe" pulse is used 
after a controllable time delay to make a measurement of 
the state of the sample at this delay. Because the delay 
between the pulses may be set very precisely, a time de- 
pendent measurement may be made on the sample. Over 
the last decade, pump-probe experiments have been used 
for the observation of Rydberg electron wavepackets [l|- 
Q , atomic spin-orbit precession Q , and molecular vibra- 
tion [a- 



Recently, the observation of spin-orbit precession in 
rubidium atoms has been demonstrated Q, and used as 
a substrate for the coherent control of the pump pulses 
1- 

In this paper, the theoretical description of spin-orbit 
precession is extended to the photoionization of long 
range Rb2 molecules. It is also shown that the same 
mechanism that enables the spin-orbit precession to be 
visible can provide a measurement of the separation of 
the atoms in the dimer. 

Two results are obtained from the extension of the 
atomic orbital precession to molecular orbital precession. 
The molecular pump-probe signal contains a high fre- 
quency component and a low frequency component, in 
contrast to the atomic case, when only a high-frequency 
component is present, at the spin-orbit precession fre- 
quency. The high frequency component decays, as a re- 
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suit of the orbital precession rate being dependent on the 
internuclear separation. The low frequency component is 
due to a vibrational oscillation of the molecule. 

Both components would be visible in a pump-probe 
ion signal. The reason is that the Hund's case (c) states 
change their admixtures of Hund's case (a) states as a 
function of internuclear separation, and the ionization 
process effectively measures the Hund's case (a) popula- 
tions rather than the Hund's case (c) populations due to 
the contrasting timescales on which the ionization pro- 
cess and spin-orbit coupling operate. This gives an in- 
ternuclear separation dependent measurement that may 
be used to observe coherent vibrational motion in the 
molecule, using a direct ionization scheme. 

The calculations of the photoionization dynamics pre- 
sented here are based on a simple model used to describe 
the rubidium atom. The model is optimized to recover 
the atomic energy levels and transition strengths of the 
rubidium atom. Ionization cross sections and the derived 
model of above threshold short pulse ionization are pre- 
sented also. 

Section [H] details the atomic model, and gives the op- 
timization process, section IIIII give the ionization cross 
sections of various states as a function of wavelength, 
section [V] gives the molecular model, and sections IVI Bl 
and IVI CI describe two experiments and their predicted 
outcomes that show the molecular orbital precession and 
the molecular vibration in the ion signal. 



II. ATOMIC MODEL 

The rubidium atom has a core containing 36 electrons 
and a single valence electron. An effective potential for 
the valence electron is employed with three parameters. 
These parameters are fit to several energy levels of the 
rubidium atom and to the transition dipolc moments of 
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several transitions. In calculating transition dipole mo- 
ments, accuracy is improved by assuming that near the 
core, the electron is effectively shielded from any incident 
radiation, and a fourth parameter is used to describe the 
shielding effect. 

The model Hamiltonian, in atomic units, is given by 

nr) = -^-{l-e-^} 2 ^. (!) 

(2) 

D(r)e(t). (3) 
(4) 



Z(r) = I + 6e~(«) C 



1 1 r)V V 2 

H(t) = V{r) + \a^l. S + \ 
D(r) = (I - e~^)r. 



Here, the first term may be considered to be the Coulomb 
potential of an effective charge, Z(r) of the core, and the 
second term is the core polarizability term. The second 
term is taken from reference [Io| , which gives the param- 
eters for the static electric core dipole polarizability, old 
and core size, p the values 8.67 E^Oq and 2.09 a respec- 
tively. The Hamiltonian contains a spin-orbit coupling 
term, a, &, c and r s are the optimised parameters of 
the model, r is the electron-nuclear separation. 1 is the 
electron orbital angular momentum operator, and s is 
the electron spin angular momentum operator, a is the 
fine-structure constant. 

The dipole operator is defined as above in order to in- 
clude a shielding effect. This can be explained physically 
by the fact that when the valence electron is in the core 
region of the atom, the core electrons will shift adiabati- 
cally in the presence of an electric field, thus shielding the 
valence electron from the electric field. It also contains 
any contribution from core rearrangement that might ac- 
company a transition in the valence electron state that 
might change the transition strength. 

The parameters given here were optimized to give a 
compromise between the energy levels and the transition 
strengths. It is possible to optimize the parameters to 
get more accurate energy levels, but at the expense of 
the transition dipole moments. 

The values for a, b and c are chosen to minimize a 
fitness function, chosen to be 



= E 



imental) ^ ' 



^] 1000 log 2 



h x THz 
Hk (model) 



Hk (experimental) 



(5) 



where the sums are over just the energy levels and tran- 
sitions selected for optimization. The fitness function 
contains a weighting between the energy level and the 
transition dipole accuracies. This value was chosen ac- 
cording to the particular requirements of the experiments 
presented here: in particular the accuracy of the energy 
levels had to be accurate to some hundreds of hxGHz. 

Values for the four parameters arc given in table U 
and the experimental observables that went into the op- 
timization are given in table [ill 



Parameter 


Value 


a/ao 


0.296760 


b 


37.68652 


c 


0.729772 


r s /a 


1.790893 



TABLE I: Optimised parameters for the atomic model pre- 
sented here. 
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E(5p 5/2 ) 


-239.454 


-239.750 


THz 


E(6p 5 /2) 


-149.939 


-150.004 


THz 


fi 5s — >5Pi/ 2 


2.99 


3.016 


q e ao 


fj, 5S — >6Pi/2 


0.236 


0.230 


q e ao 


fi 5S^7P 1/2 


0.0813 


0.0762 


qedo 


H 5S^8P 1/2 


0.0407 


0.0382 


q e ao 


A* 5S^9P 1/2 


0.0252 


0.0234 


q e ao 


fi 5S — s-5P 3/2 


2.98 


3.008 


q e ao 


fi 5S — s-6P 3/2 


0.255 


0.265 


q e ao 


fi 5S — ^7P 3/2 


0.0950 


0.0963 


q e ao 


fi 5S^8P 3/2 


0.0504 


0.0517 


q e ao 


fi 5S^9P 3/2 


0.0326 


0.0333 


q e ao 


fi 5P 3/2 — >6d 5/2 


0.677 


0.631 


q e ao 


AEso 


7.123 


7.513 


THz 



TABLE II: The experimental observables that went into the 
optimizations. The energies of 18 levels were optimized, as 
were eleven transition dipole moments. 



The model is much simpler than multi-electron mod- 
els which require configuration interaction calculations or 
self consistent methods to be applied, and allow calcula- 
tions to be performed much more easily. When compared 
to other single-electron models, such as reference [l2[, the 
model presented here recovers remarkably accurate tran- 
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FIG. 2: The continuous-wave photoionization cross section as 
a function of wavelength for the ground state rubidium atom. 
The blue line is the model's predictions. The point (a) shows 
the result from reference [ljj], including error. The two re- 
sults are i n ag reement. The points (b) shows the results from 
reference, [lj]. These are a relative magnitude measurement 
and so have been scaled appropriately. 



FIG. 1: (Color online) The energy levels of the rubidium 
atom. The atoms in the experiments discussed here are ini- 
tially in the ground state, 5s. A short pulse transfers them 
to 5p states. An ionizing pulse then transfers population to 
continuum states with either s or d symmetry. If both acces- 
sible 5p states are populated and the ionizing pulse is short 
enough, then interference effects cause the precession of the 
orbital angular momentum about the total angular momen- 
tum to be visible in the ion signal. 

sition strengths and ionization cross sections, as well as 
treating spin-orbit coupling accurately. These are essen- 
tial for estimating the visibility of the spin-orbit preces- 
sion in an ion signal, and so the model is well suited to 
the applications presented here. 

The numerical approach to solving the Schrodinger 
equation for Eqn. [3] is given in the Appendix. 



III. ATOMIC PHOTOIONIZATION CROSS 
SECTION 

Using this atomic model, the photon energy dependent 
photoionization cross section of the ground state is cal- 
culated for verification as described in the Appendix. 

The continuous wave photoionization cross section as a 
function of photon energy for the ground state is shown 
in Fig. [2] The model agrees with experiment [ill 13 
for the photoionization of the ground state, but overesti- 
mates cross sections for the photoionization cross section 
of the 5P state by around 20% when compared to experi- 



ment [15(. This could be for a variety of reasons. Firstly, 
the error bars on the experimental result are wide, and 
it is possible that the true value is quite close to the 
calculation presented here. Secondly, the photionization 
cross sections of the 5p states depend on the transition 
dipolc moments from p orbitals to d orbitals. Due to the 
scarcity of experimental transition dipolc moments from 
p to d orbitals, the model was only optimized to recover 
a single p to d transition. Therefore it is quite likely that 
although the s and hp electron wavefunctions are quite 
well recovered, the d orbital wavefunctions could be less 
accurate, and this affects the p orbital ionization cross 
sections. 

For short-pulse photoionization, the ionization prob- 
ability may be expressed as the magnitude of the final 
wavefunction when projected onto the scattering elec- 
tronic states: 

P(ion) ~ |p ion £/Vi| 2 , (6) 

where U is the propagator from before the start of the 
pulse to after the end of the pulse, pj on projects onto 
only ionised electronic states, and ip[ is the initial state. 
Equivalcntly, this can be expressed as the expectation 
value of an operator, M . 

The atom has six 5P excited states: two spin states 
multiplied by the three orbital angular momentum pro- 
jection quantum numbers: -1, and 1. The axis of pro- 
jection will be chosen to be the electric field axis of the 
ionizing electric field. In the fine structure basis, the 
5Pi/2 state has 2-fold degeneracy and the 5P 3 / 2 state has 
4-fold degeneracy. Any operator acting on this basis may 
be expressed as a six by six matrix, but the reversal sym- 
metry ( mi — > —mi, m s — > —m s ) means that the states 
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with positive total angular momentum are never coupled 
to states with negative angular momentum, and so we 
can limit the discussion to the positive to positive ma- 
trix elements in the understanding that the negative to 
negative elements are the same. More generally, neither 
the electric field nor the spin-orbit coupling changes the 
total projected angular momentum of the atom, so the 
four possible values, -3/2, -1/2, +1/2 and +3/2 remain 
separate: Nothing in this model couple between values 
for this quantum number. The ionization operator will 
therefore be block diagonal with each block representing 
a different total angular momentum. 

The three positive total angular momentum basis 
states are 



\J- 
\J- 
\J 



1/2, m, = 1/2), 
3/2, mj = 1/2), 
= 3/2, mj = 3/2). 



(7) 

(8) 
(9) 



If an arbitrary 5P wavefunction \tp) is represented by a 
vector a: 

\i>) = ai |i = l/2,mi = 1/2) 
+a 2 \j = 3/2, mj = 1/2) 
+a 3 |j = 3/2,m i = 3/2), 

then the ionization operator, A/, whose expectation value 
gives the ionization probability is a three by three Hcr- 
mitian matrix. 

As an example, the ionization operator whose expecta- 
tion value gives the ionization probability for a Gaussian 
probe pulse, centred at 650 THz, with a full width at half 
maximum of 18 THz, and with a total fluence of 0.3 Jm -2 
is given as 



M 



( 0.997 -0.223 



V 




(10) 



A basis change can be made to the s,m s ,l, rri£ basis, 
so that the 5P state is expressed as 

\i>) = b 1 \m e = 0,m s = +l/2) 
+b 2 \m e = l,m s = -1/2) 
+b 3 \m e = l,m s = +1/2). 

In this basis, the ionization matrix takes the form 



M 



' 1.355 0.0297 

0.0297 0.860 I x 10~ 2 
v 0.817 



(11) 



The small off-diagonal elements show that the ioniza- 
tion process effectively measures the populations in the 
mi basis. The physical reason for this is that the pulse 
is shorter than the spin-orbit interaction time, and so 
the electron spins may be neglected, leaving the differ- 
ing £, rri£ quantum numbers as the elements that affect 



the ionization probability. The difference in ionization 
probability between the two mg = 1 states is due to the 
subtly different spatial wavefunction of the j = 3/2 and 
j = 1/2 states. 

The ionization matrix in this example, and in all the 
cases presented, was calculated by propagating the three 
initial states under the influence of the electric field and 
using their ionization probability and the inner product 
between the resulting ionized wavefunctions to infer ma- 
trix elements. Although the field used here is in the weak 
field regime, this method of parameterizing the effect of 
the field applies to any field. 

The ionization process will be approximated in the 
molecular simulations as 



M 



\rw = 0)1.36- 10" 



(me 



0| 



\m e = 1)0.84- l0~ A (m e = 1| 



(12) 
(13) 
(14) 



The approximation is correct to within a few percent, 
which is also roughly the error in transition dipole mo- 
ments of the atomic model. 



A. Ionization operator dependence 

The framework described above is used to find the 
ionization parameters as a function of pulse bandwidth 
and centre frequency. The two quantities that determine 
whether the ionization operator is selective for one orbital 
orientation are the angle between the operator's principle 
axis and the \mt = 0) axis, and the ratio of ionization 
probabilities for the \mg = 1) and \rri£ = 0) states. These 
two values are shown in Fig. |3] as a function of probe 
wavelength and bandwidth. The ionization operator can 
be seen to vary slowly with laser parameters over these 
intervals, becoming more orbital-alignment selective at 
smaller probe bandwidths. 



IV. ATOMIC SPIN-ORBIT DYNAMICS 

The spin-orbit dynamics of the first excited state of the 
rubidium atom are of most interest here due to the ease 
of production of the state. The two levels in question are 
illustrated in Fig. [1] 

The atom has a single valence electron, giving it an 
electron spin of i. In the first excited state, the 5P state, 
the electron orbital has one unit of angular momentum. 
This gives the atom 6 allowed angular momentum states. 
Two pairs of these are coupled by the spin-orbit coupling: 



\mi 



0,m s = ±i) is coupled to \mi = ±1, 



Since the atom is transferred by a short laser pulse to 
the \mi = 0) state, where the quantization axis is chosen 
along the electric field of the exciting laser pulse, the spin- 
orbit coupling causes an oscillation between the \mt = 0) 
and \mi = 1) states. Since the \me = state is about 
60% easier to ionize than the \m,£ = 1) state, this gives 
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FIG. 3: The ionization operator is a function of the probe 
pulse parameters. The operator has a principle axis defined 
as the state that is most easily ionized. The angle between 
this and the \mi = 0} state is shown in the first panel. The 
ratio of ionization probabilities for the \mi = 1} compared to 
the \mt = 0) state is shown in the second panel. The result 
is that both quantities have a weak dependence on the probe 
pulse, but narrower bandwidth probes are more sensitive to 
orbital alignment. 



5s 5p (M^-=0) 5p(M,-=l) 




FIG. 5: The rubidium ground state and first excited states. In 
the experiments discussed, the atom starts off in the ground 
state (5s) and is transferred to the first excited state. If the 
transition is instantaneous, then the excited population will 
be in the mj = state (centre). This is not an eigenstate 
of the full Hamiltonian since it is spin orbit coupled to the 
mi = 1 state (right). The result is that the electron orbital 
oscillates between the me = and mi = 1 states. Since 
the two states have different ionization probabilities when an 
ionizing pulse is applied, the probability of ionizing the atom 
oscillates. The surfaces are a contour of p(r)|r| 2 . 



the pump-probe signal of the atom an oscillation at the 
spin-orbit frequency of 7.123 THz. 

The initial state and the \mg = 0) and \mp = 1) states 
are shown in Fig. O 




FIG. 4: The wavefunction of the valence electron at the peak 
of a short ionizing laser pulse. Contours of probability density 
multiplied by the electron-core separation squared are chosen 
at intervals to highlight the exiting electron's wavefunction, 
which is a superposition of s and d waves. Because of the 
short timescale, the electron has not had time to travel far 
from the core. The shelled structure is due to the electron 
being emitted at a specific phase of the incident radiation. 
The electric field axis is aligned vertically. 



V. EXTENSION TO MOLECULAR CASE 

The calculations can be extended to the molecular case 
by increasing the 8 atomic states (2 ground, 6 excited) to 
15 molecular states (3 ground, 12 excited), which make 
up all the accessible electronic states of the molecule, 
assuming it only absorbs a single photon and is initially 
in a triplet state. The intcrnuclcar separation is added 
as a continuous variable, and the angle the internuclear 
axis makes with the electric field of the laser pulses is 
added as a constant parameter since it does not change 
appreciably over the timescale of the experiment. The 
Schrodinger equation is thus converted from a discrete 
8 by 8 matrix equation into a one dimensional partial 
differential equation with 15 channels. The details of the 
model are given in a previous publication [i"6j . 

The ionization step is modelled by projecting the 
molecular state onto the separated atomic states and 
taking the expectation value of the ionization matrix 
(Eqn. fl~3|) as the ionization probability. 

VI. EXPERIMENTS 

The application of the atomic model to the study of the 
direct photoionization of long range molecules is demon- 
strated by considering two experiments. 

The two experiments outlined here arc both possible, 
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and an effort has been made to ensure they are feasi- 
ble with current technology. Both have the same initial 
state which can be prepared using a ramped magnetic 
field over a Feshbach resonance. In both experiments, a 
pump pulse is applied, then after a delay, a probe pulse 
with the same polarization is applied. The probe pulses 
are identical, but the pump pulses arc different. The po- 
larization was chosen parallel to the magnetic field used 
to associate the molecules. 

In the first experiment, the pump pulse is a Fourier 
limited Gaussian pulse whose bandwidth spans both the 
transition from the ground state, 5S, to the 5P]/2 and 
to the 5P3/2 state. This maximises the visibility of the 
spin-orbit precession. In the second experiment, the same 
Gaussian pump pulse is spectrally cut to remove inten- 
sity at the 5S — 5P 3 / 2 transition. This prevents the 
preparation of a precessing orbital. As a result, the only 
mechanism by which the orbital can change orientation is 
through the interaction of the two atoms in the molecule. 
Since this depends on internuclear separation, the oscil- 
lation in the internuclear separation coordinate becomes 
apparent in the pump-probe signal. 



A. Initial state 

The initial state for each experiment is prepared using 
a magnetic field sweep across a Feshbach resonance in a 
cold 87 Rb gas with a high phase space density in order 
to associate atom pairs to loosely bound molecules. 

The ground electronic state rubidium dimer has a dis- 
crete number of bound states plus several continua of 
scattering states, each with a different electron and nu- 
clear spin state. The energies of the bound and contin- 
uum states change as a function of magnetic field and 
spin state. The dynamics and molecular properties of 
long range dimers bound by less than a few hxGHz is 
dictated by the hyperfine and Zeeman interactions. Sev- 
eral reviews [l7j | and other publications are concerned 
with the behaviour of such molecules, which constitute a 
major area of study in their own right. 

The use of a Feshbach resonance to create molecules 
was demonstrated in 2002 [lj|. Work had already been 
done on using a magnetic field to control atom-atom in- 
teractions both theoretically and experimentally [l9|, [2(| ■ 
This mechanism is the mechanism for creating the initial 
state in the work here. An essay by D. Kleppner presents 
a fuller picture of the use of Feshbach resonances in this 
context 2l|. Magnetic sweep association of the species 



under the conditions studied here has been demonstrated 
experimentally (22j . 

Two interacting atoms are prepared in their high-field 
seeking / = 1 state. The spin-spin interaction between 
the two atoms is neglected, which constrains the atom 
pair to stay in the mj = 2 state, since the remaining hy- 
perfine, Zeeman, and electronic interactions do not break 
the rotational symmetry of the electron plus nuclear spins 
about the magnetic field axis. For numerical reasons, the 
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FIG. 6: The highest bound state energies of the rubidium 
dimer as a function of magnetic field. For the experiments 
outlined here, the atom pairs are initially unbound at a mag- 
netic field of greater than 1007 G (a) . The magnetic field is 
swept down to 940 G (b) for the first experiment, and to 850 G 
(c) for the second experiment. This prepares a suitable initial 
state for each experiment. 



internuclear separation is constrained to be in a box. The 
highest 9 vibrational states and lowest 191 box states for 
each of the 5 accessible spin states are taken as a ba- 
sis set. The Hamiltonian, consisting of the vibrational 
energies, the electron and nuclear Zeeman energies, and 
the hyperfine interaction, was diagonaliscd for each mag- 
netic field. The energy levels of the dimer as a function 
of magnetic field for this model are shown in Fig. [5] 

By following a single state adiabatically from above 
1007 G, where it is a box state, to 940 G, where it is bound 
by 24 hxMHz, or to 850 G for the second experiment, a 
suitable initial state is prepared. 

For the 940 G state, the state is well represented by the 
highest single channel vibrational state (whether singlet 
or triplet) multiplied by a spin state l^) which is 74% 
\S = I, M s = -1,1 = 3, Mj = 3), 

12% \S = 1,M 3 = 0,1 = 3, Mi = 2), and 12% 
\S = 0,M S = 0,1 = 2, Mi = 2). The superposition 
is coherent, but the phases between these three states 
are not measured by the experiment, since the rest of 
the experiments take place on a timescale too short to 
change the nuclear spins and so the nuclear spins effec- 
tively measure the initial state of the electron spin. For 
this state, the singlet population is neglected. There are 
two justifications for this. The first is that it is only 12% 
of the population, and the second is that it is expected 
to behave similarly to the triplet population: It will spin- 
orbit precess at first, but molecules at different internu- 
clear separations will spin-orbit precess at different rates, 
so the visibility of the precession will wash out. 

For the 850 G state, the state is well represented by 
the fifth highest triplet state with 70% in the 
\S = 1,M S = 1,1 = 3, Mi = 1) state, and 28% in the 
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FIG. 7: The ionization probability of a long-range molecule 
prepared as described in section IVI Al at 940 G in a pump- 
probe experiment where the molecule is excited to the states 
associated with the 5S + 5P asymptote and then ionized 
directly. The pump and probe are sub-picosecond pulses. 
This is the calculated experimental outcome for experiment 1, 
where the molecules are put in a superposition of fine struc- 
ture states in order to cause spin-orbit precession. The spin- 
orbit precession is clearly visible in the signal here. 



\S = 1, M s = -1,1 = 3, Mi = 3) state. 



B. Experiment 1: Broadband uncut pump 



FIG. 8: The relative number of molecular ions yielded as a 
function of internuclear separation for a long-range molecule 
prepared as described in section [VI Al in a pump- probe exper- 
iment where the molecule is excited to the states associated 
with the 5S + 5P asymptote and then ionized directly. The 
angle between the internuclear axis and the electric field is 1.4 
radians. The contribution is shown as a function of internu- 
clear separation and time. This is the population density mul- 
tiplied by ionization probability for experiment 1. It shows the 
spin orbit precession, and how the contributions from differ- 
ent internuclear separations become dephased leading to the 
signal washing out. The figure has been deliberately under- 
sampled to increase the apparent oscillation period to around 
lps rather than 1/7 ps. 



In this experiment, we wish to observe spin-orbit pre- 
cession in the long range Rb2 dimer. The molecule is 
prepared as described in section IVl Al This initial state 
is very loosely bound, which increases the dephasing time 
for the spin-orbit precession, allowing a clearer signal. 
The molecules are excited with a Gaussian pump probe 
with a centre wavelength of 375 THz and a full width 
at half maximum intensity of 30 THz. After a variable 
time delay, the molecules are ionized by a second Gaus- 
sian laser pulse, centred at 650 THz, with a full width 
at half maximum of 18 THz, and with a total fluence of 
0.3 Jm~ 2 . The ionization probability as a function of de- 
lay is measured by counting the number of molecular ions 
produced. 

The experiment is modelled by taking the initial state, 
propagating it under the influence of the first pulse. At 
each timestep, the expectation of the ionization operator 
can be taken, which gives the ionization probability of 
the molecule at the specified delay. The calculation is 
repeated for 21 different angles between the internuclear 
axis and the electric field and averaged. This ionization 
probability is given in Fig. [7] The contribution to this 
signal as a function of internuclear separations for one 
angle is given in Fig. [8] 

The spin-orbit precession is initially visible as a 7 THz 
oscillation, with a good signal to background, but then 



decays due to the interaction of the two atoms. 



C. Experiment 2: Spectrally cut pump 

In this experiment, the aim is to observe the internu- 
clear oscillations rather than spin-orbit precession. As 
a result, a different initial state is chosen that has the 
fifth highest triplet state occupied rather than the high- 
est. This shifts the molecular population to separations 
where the ionization process is more sensitive to separa- 
tion, and where the population oscillates more uniformly. 
The pump pulse is spectrally cut to suppress the orbital 
precession. The probe is the same as for experiment 1. 
The experiment is modelled as before, and the expected 
ionization probability is given in Fig. |H1 The contribu- 
tions to this signal are shown in Fig. [lOj 

The spin-orbit precession is suppressed as hoped. A 
persistant long period oscillation may be observed mixed 
with a short period oscillation. The long period signal is 
due to the oscillation in the internuclear separation. The 
short period signal is due to spin-orbit precession. 



8 



x 10" 




150 



Time (ps) 



FIG. 9: The ionization probability of a long-range molecule 
prepared as described in section IVI Al at 850 G in a pump- 
probe experiment where the molecule is excited to the states 
associated with the 5S + 5P asymptote and then ionized di- 
rectly. This is the signal for the second experiment. The 
pump pulse is spectrally cut to suppress the spin-orbit pre- 
cession. As a result, and as a result of the different initial 
state, the molecular dynamics are visible as a long period 
contribution to the signal. 
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tion in the ion signal washes out depends on the vibra- 
tional state or states that the molecules are in. 

The second experiment is similar to the first exper- 
iment discussed in reference (l6| . In the earlier work, 
it was shown that a coherently oscillating excited state 
wavefunction could be formed in long-range molecules. 
It was stated that if a position sensitive measurement 
could be made, then the molecular dynamics could be 
observed. In this paper, the result is extended by show- 
ing that a position dependent measurement can indeed 
be made, using orbital alignment. The initial state has 
to be modified from the earlier work, since it is difficult 
experimentally to use an above threshold ionization laser 
on a gas in a magneto-optical trap due to atomic ioniza- 
tion. In this way, the pump-probe experiment can now 
be modelled from initial state preparation to the creation 
of the molecular ions. The observation of wavepacket dy- 
namics in the excited state of long range molecules would 
be a step towards the control and stabilization of such 
molecules, and would be a valuable contribution to the 
field of ultracold chemistry. 



VII. CONCLUSIONS 

A simple atomic model is presented that is accurate 
enough to model spin-orbit precession and photoioniza- 
tion. A method for extending atomic spin-orbit preces- 
sion calculations to the case of long range molecules is 
presented. This is demonstrated by the calculation of 
the ionization signal for two pump-probe experiments. 
A consequence of this is that it is shown that long range 
molecules can exhibit spin-orbit precession, but that it 
washes out. A second consequence is that it is shown that 
the mechanism that causes spin-orbit precession to affect 
the ionization signal also causes an internuclear separa- 
tion oscillation to affect the ionization signal, allowing 
pump-probe experiments to observe coherent wavepacket 
behaviour in ultracold long range molecules. 
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FIG. fO: The contribution to the total ion yield of a long- 
range molecule prepared as described in section IVI Al as for 
Fig. but for experiment 2. It shows the internuclear sep- 
aration oscillations clearly. The admixture of excited states 
that ionize easily change with internuclear separation, and 
this results in a time-dependent ion signal. 
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APPENDIX 



D. Discussion 

The first experiment would demonstrate that it is 
possible to observe spin-orbit precession in long-range 
molecules. One application for this is to verify which 
molecular vibrational state has been prepared by a mag- 
netic field sweep, since the speed with which the oscilla- 



The eigenstates and eigenenergies of the Schrddingcr 
equation are found in the basis defined by the quantum 
numbers of total angular momentum j = I + s, the pro- 
jected total angular momentum mj = mi + m s , and the 
electron orbital angular momentum I — a basis for which 
equation [3] is diagonal. The lowest fOO energy levels for 
the s,p,d,f and g orbitals are used. The calculations are 
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performed in a spherical box of 1600 ciq. This ensures 
that the lowest 100 energy levels contain scattering states 
with an energy range that covers the possible free elec- 
tron energies in the experiments discussed here. The fi- 
nite spacing of the scattering states puts an upper bound 
on the duration of any ionizing pulse that can be accu- 
rately modelled. However, all the pulses used here are 
short enough that the scattering states included in the 
basis set provide an accurate representation of the con- 
tinuum. The energy levels were obtained using fourth 
order Runge-Kutta. 

The time-dependent Schrodinger equation was solved 
using a matrix representation of the propagator for each 
timestep, chosen to be shorter than the time variation 
in the Hamiltonian. For weak fields, a 1st or 2nd order 
pcrturbative expansion may be used. 



Continuous- wave photoionization cross sections can be 
calculated in this framework using the equation 



where /ij is the transition dipolc moment from the initial 
state to the final scattering state at energy E normalized 
to a box of length L, and u is the angular frequency of 
the incident linearly polarized light. The sum is over the 
various electronic configurations of the scattering state: 
the different allowed orbital and spin angular momenta 
of the exiting electron. The equation becomes exact in 
the limit that L — > oo. 
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